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Abstract 
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■ We study polar ized-spin structure functions of the nucleon within the 

£S| ■ bosonized Nambu-Jona-Lasinio model where the nucleon emerges as a chi- 

C\ . 

ral soliton. We present the electromagnetic polarized structure functions, 



gi (x) and g<i (x) for ep scattering and discuss various sum rules in the va- 
lence quark approximation. This approximation is justified because in this 

a.: 

model axial properties of the nucleon are dominated by their valence quark 



X 



contributions. We find that these structure functions are well localized in 
the interval < x < 1. We compare the model predictions on the polarized 
structure functions with data from the E143 experiment by evolving them 
from the scale characteristic of the NJL-model to the scale of the data. Addi- 
tionally a comparison is made with parameterized data at a momentum scale 
commensurate with the model calculation. 
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I. INTRODUCTION 



Over the past decade, beginning with the measurement of nucleon spin-polarized struc- 
ture function, gi(x, Q 2 ) by the EMC Jl|] at CERN and most recently with the spin-structure 
function g 2 (x,Q 2 ) in the E143 experiment |§ at SLAC, a wealth of information has been 
gathered on the spin-polarized structure functions of the nucleon and their corresponding 
sum rules (see in addition ||, ||, 0, ||, [|7|], ||). Initially the analysis of these experiments 
cast doubt on the non-relativistic quark model || interpretations regarding the spin content 
of the proton. By now it is firmly established that the quark helicity of the nucleon is much 
smaller than the predictions of that model, however, many questions remain to be addressed 
concerning the spin structure. As a result there have been numerous investigations within 
models for the nucleon in an effort to determine the manner in which the nucleon spin is 
distributed among its constituents. One option is to study the axial current matrix elements 
of the nucleon such as (N\A^\N) = 2Aq i S fl , which, for example, provide information on the 
nucleon axial singlet charge 

g\ = (N\A° 3 \N) = (An + Ad + As) = T^Q 2 ) + r™(Q 2 ) . (1) 

Here Aq are the axial charges of the quark constituents and T^(Q 2 ) = dxg^{x, Q 2 ) is 
the first moment of the longitudinal nucleon spin structure function, g^ (x, Q 2 ). Of course, 
it is more illuminating to directly compute the longitudinal and transverse nucleon spin- 
structure functions, g\ (x, Q 2 ) and gr(x, Q 2 ) = gi(x, Q 2 )+g 2 (x, Q 2 ), respectively as functions 
of the Bjorken variable x. We will calculate these structure functions within the Nambu- 



Jona-Lasinio (NJL) [JTOj] chiral soliton model [11 



Chiral soliton models are unique both in being the first effective models of hadronic 
physics to shed light on the so called "proton-spin crisis" by predicting a singlet combination 



in accord with the data [12], and in predicting a non-trivial strange quark content to the 
axial vector current of the nucleon fll2| , fI5|| , [[14]], fI5fl ; about 10 — 30% of the down quarks 
(see |16| and [17] for reviews). However, while the leading moments of these structure 



functions have been calculated within chiral soliton models, from the Skyrme model |TH 



TT| and its various vector-meson extensions, to models containing explicit quark degrees of 



freedom such as the (NJL) model JIU], the nucleon spin-structure functions themselves have 



not been investigated in these models. Soliton model calculations of structure functions 



were, however, performed in Friedberg-Lee |HJ and color-dielectric 0] models. In addition, 
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structure functions have extensively been studied within the framework of effective quark 
models such as the bag-model |£| , and the Center of Mass bag model |23| . These models 
are confining by construction but they neither contain non-perturbative pseudoscalar fields 
nor are they chirally symmetrical. To this date it is fair to say that many of the successes of 
low-energy effective models rely on the incorporation of chiral symmetry and its spontaneous 
symmetry breaking (see for e.g. ]j26|| ). In this article we therefore present our calculation 
of the polarized spin structure functions in the NJL chiral soliton model [27]], [26||. Since 
in particular the static axial properties of the nucleon are dominated by the valence quark 
contribution in this model it is legitimate to focus on the valence quarks in this model. 

At the outset it is important to note that a major difference between the chiral soliton 
models and models previously employed to calculate structure functions is the form of the 
nucleon wave-function. In the latter the nucleon wave-function is a product of Dirac spinors 
while in the former the nucleon appears as a collectively excited (topologically) non-trivial 
meson configuration. 

As in the original bag model study ]^9[] of structure functions for localized field configura- 
tions, the structure functions are most easily accessible when the current operator is at most 
quadratic in the fundamental fields and the propagation of the interpolating field can be 
regarded as free. Although the latter approximation is well justified in the Bjorken limit the 
former condition is difficult to satisfy in soliton models where mesons are fundamental fields 
(e.g. the Skyrme model |TB[, |TP[, the chiral quark model of ref. or the chiral bag model 
|3l|). Such model Lagrangians typically possess all orders of the fundamental pion field. In 
that case the current operator is not confined to quadratic order and the calculation of the 
hadronic tensor (see eq. (@) below) requires drastic approximations. In this respect the chi- 
rally invariant NJL model is preferred because it is entirely defined in terms of quark degrees 
of freedom and formally the current possesses the structure as in a non-interacting model. 
This makes the evaluation of the hadronic tensor feasible. Nevertheless after bosonization 
the hadronic currents are uniquely defined functionals of the solitonic meson fields. 

The paper is organized as follows: In section 2 we give a brief discussion of the standard 
operator product expansion (OPE) analysis to establish the connection between the effective 



1 In the cloudy bag model the contribution of the pions to structure functions has at most been 



treated perturbatively [24], [25|. 
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models for the baryons at low energies and the quark-parton model description. In section 
3 we briefly review the NJL chiral soliton. In section 4 we extract the polarized structure 
functions from the hadronic tensor, eq. ([16]) exploiting the "valence quark approximation" . 
Section 5 displays the results of the spin-polarized structure functions calculated in the 
NJL chiral soliton model within this approximation and compare this result with a recent 
low-renormalization point parametrization JSJ. In section 6 we use Jaffe's prescription 
|33| to impose proper support for the structure function within the interval x G [0,1]. 
Subsequently the structure functions are evolved ||34|| , fl35| , [j36fl from the scale characterizing 
the NJL-model to the scale associated with the experimental data. Section 7 serves to 
summarize these studies and to propose further explorations. In appendix A we list explicit 
analytic expressions for the isoscalar and isovector polarized structure functions. Appendix 
B summarizes details on the evolution of the twist-3 structure function, g 2 (x, Q 2 ). 



II. DIS AND THE CHIRAL SOLITON 



It has been a long standing effort to establish the connection between the chiral soliton 
picture of the baryon, which essentially views baryons as mesonic lumps and the quark parton 
model which regards baryons as composites of almost non-interacting, point-like quarks. 
While the former has been quite successful in describing static properties of the nucleon, 
the latter, being firmly established within the context of deep inelastic scattering (DIS), has 
been employed extensively to calculate the short distance or perturbative processes within 
QCD. In fact this connection can be made through the OPE. 

The discussion begins with the hadronic tensor for electron-nucleon scattering, 

wm = ^J*t j q 'Hm 4(o)} \n) , (2) 

where J M = ?(^)7 M Q?(^) is the electromagnetic current, Q = ( |, is the (two flavor) quark 
charge matrix and \N) refers to the nucleon state. In the DIS regime the OPE enables one to 
express the product of these currents in terms of the forward Compton scattering amplitude 
T l _ iv (q) of a virtual photon from a nucleon 

T, u (q) = e«*(N\T (J„(0 Jj(0)) \N) , (3) 
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by an expansion on the light cone (£ 2 — > 0) using a set of renormalized local operators [37 



In the Bjorken limit the influence of these operators is determined by the twist, r or the 
light cone singularity of their coefficient functions. Effectively this becomes a power series 
in the inverse of the Bjorken variable x = —q 2 /2P-q, with P^ being the nucleon momentum: 

= E C-X eiu(l,P,S) CUQV^,a s (^))O n T ^ 2 )(^ 1 ■ (4) 

n,i,r VX/ ^ 

Here the index % runs over all scalar matrix elements, 0™j(/i 2 ), with the same Lorentz 
structure (characterized by the tensor, e 1 ). Furthermore, S 11 is the spin of the nucleon, 
(S 2 = —1 , S ■ P =0) and Q 2 = —q 2 > 0. As is evident, higher twist contributions are 
suppressed by powers of l/Q 2 - The coefficient functions, C^^Q 2 / fi 2 ,a s (fi 2 )) are target 
independent and in principle include all QCD radiative corrections. Their Q 2 variation is 
determined from the solution of the renormalization group equations and logarithmically 
diminishes at large Q 2 . On the other hand the reduced-matrix elements, 0™j(/z 2 ), depend 
only on the renormalization scale /i 2 and reflect the non-perturbative properties of the 
nucleon |39] . 



The optical theorem states that the hadronic tensor is given in terms of the imaginary 
part of the virtual Compton scattering amplitude, W^ u = t^Thi T^ v . From the analytic 
properties of T fiu (q), together with eq. (^) an infinite set of sum rules result for the form 
factors, Wi(x, Q 2 ), which are defined via the Lorentz covariant decomposition W fiu (q) = 
e 1 Wi (x, Q 2 ). These sum rules read 

£dx x- 1 (y,Q 2 ) = £ CI, (Q 2 /^ 2 ,^ 2 )) 0^){^ - 1 . (5) 

In the impulse approximation (i.e. neglecting radiative corrections) [[E|-|42| one can directly 
sum the OPE gaining direct access to the structure functions in terms of the reduced matrix 
elements C"j(/^ 2 ). 

When calculating the renormalization-scale dependent matrix elements, C^j(/i 2 ) within 
QCD, fi 2 is an arbitrary parameter adjusted to ensure rapid convergence of the perturbation 
series. However, given the difficulties of obtaining a satisfactory description of the nucleon 
as a bound-state in the Q 2 regime of DIS processes it is customary to calculate these matrix 
elements in models at a low scale fi 2 and subsequently evolve these results to the relevant 
DIS momentum region of the data employing, for example, the Altarelli-Parisi evolution 



34j, [pq| . In this context, the scale, /r ~ Aq CD , characterizes the non-perturbative regime 
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where it is possible to formulate a nucleon wave-function from which structure functions are 
computed. 

Here we will utilize the NJL chiral-soliton model to calculate the spin-polarized nucleon 
structure functions at the scale, /i 2 , subsequently evolving the structure functions according 
to the Altarelli-Parisi scheme. This establishes the connection between chiral soliton and 
the parton models. In addition we compare the structure functions calculated in the NJL 
model to a parameterization of spin structure function [|32| at a scale commensurate with 
our model. 



III. THE NUCLEON STATE IN THE NJL MODEL 



The Lagrangian of the NJL model reads 

£ = q(i$ - m°)q + 2G NJL £ [[q T -qf + (g^ 7 5<?) 2 ) . (6) 

Here q, rh° and Gnjl denote the quark field, the current quark mass and a dimension- 
ful coupling constant, respectively. When integration out the gluon fields from QCD a 
current-current interaction remains, which is meditated by the gluon propagator. Replac- 
ing this gluon propagator by a local contact interaction and performing the appropriate 
Fierz-transformations yields the Lagrangian (||) in leading order of 1/N C |43|| , where N c 
refers to the number of color degrees of freedom. It is hence apparent that the interaction 
term in eq. (H) is a remnant of the gluon fields. Hence gluonic effects are included in the 
model described by the Lagrangian @. 

Application of functional bosonization techniques |44| to the Lagrangian (||) yields the 
mesonic action 



A = Tr A log(iD) + — - — f d 4 x tr (m° (M + M f ) - MM f ) , (7) 
D = 0-(M + M f ) - 7 5 (M - M f ) . (8) 



The composite scalar (S) and pseudoscalar (P) meson fields are contained in M = S+iP and 
appear as quark-antiquark bound states. The NJL model embodies the approximate chiral 
symmetry of QCD and has to be understood as an effective (non-renormalizable) theory of 
the low-energy quark flavor dynamics. For regularization, which is indicated by the cut-off 
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A, we will adopt the proper-time scheme [45] . The free parameters of the model are the 



current quark mass m°, the coupling constant Gnjl and the cut-off A. Upon expanding 
A to quadratic order in M these parameters are related to the pion mass, = 135MeV 
and pion decay constant, f n = 93MeV. This leaves one undetermined parameter which we 
choose to be the vacuum expectation value m = (M). For apparent reasons m is called 
the constituent quark mass. It is related to m°, Gnjl and A via the gap-equation, i.e. the 



equation of motion for the scalar field S ||44|| . The occurrence of this vacuum expectation 



value reflects the spontaneous breaking of chiral symmetry and causes the pseudoscalar fields 
to emerge as (would-be) Goldstone bosons. 



As the NJL model soliton has exhaustively been discussed in recent review articles |26 
46fl we only present those features, which are relevant for the computation of the structure 
functions in the valence quark approximation. 

The chiral soliton is given by the hedgehog configuration of the meson fields 

Mh(cc) = m exp (it ■ xQ(r)) . (9) 

In order to compute the functional trace in eq. (^) for this static configuration we express 
the Dirac operator (§) as, D = 270 (<9* — h) where 

h = ex ■ p + m exp (275T • xQ(r)) (10) 

is the corresponding Dirac Hamiltonian. We denote the eigenvalues and eigenf unctions of 
h by e M and ty^, respectively. Explicit expressions for these wave-functions are displayed 
in appendix A. In the proper time regularization scheme the energy functional of the NJL 



model is found to be 27,26 



£[e] = ^ v (l + sgn(e v )) +^ r^£=5>p(-^) 
2 2 J 1 a 2 \/47rs 3 V v ' 



/A 2 VATCS C 

+ m 2 JlJ tfV(l-cos6(r)), (11) 



with Nc = 3 being the number of color degrees of freedom. The subscript "v" denotes the 
valence quark level. This state is the distinct level bound in the soliton background, i.e. 
—m < e v < in. The chiral angle, 0(r), is obtained by self-consistently extremizing E[Q] 
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States possessing good spin and isospin quantum numbers are generated by rotating the 



hedgehog field |T9 
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Af (as, *) = A(t)M K (x)A*(t) , 



(12) 



which introduces the collective coordinates A(t) e SU(2). The action functional is expanded 



1 27 1 in the angular velocities 

2A ] (t)A(t) =ir-Q . (13) 
In particular the valence quark wave-function receives a first order perturbation 

v& v (*, t) = e-^A(t) ( * v (aj) + W =: e^^)^)- (14) 

Here ip Y {x) refers to the spatial part of the body-fixed valence quark wave-function with the 
rotational corrections included. Nucleon states \N) are obtained by canonical quantization 
of the collective coordinates, A(t). By construction these states live in the Hilbert space of 
a rigid rotator. The eigenf unctions are Wigner D-functions 

(A\N) = ±D%_ Ja (A) , (15) 

with 7 3 and J3 being respectively the isospin and spin projection quantum numbers of the 
nucleon. 



IV. POLARIZED STRUCTURE FUNCTIONS IN THE NJL MODEL 



The starting point for computing nucleon structure functions is the hadronic tensor, 
eq. @. The polarized structure functions are extracted from its antisymmetric piece, 
Wffl = (W IJU , — W vli )/2i. Lorentz invariance implies that the antisymmetric portion, charac- 
terizing polarized lepton-nucleon scattering, can be decomposed into the polarized structure 
functions, gi(x,Q 2 ) and g 2 (x,Q 2 ), 

W$(q) = ie^^f j#iO, Q 2 )S° + (s° - ^-V) <fc(s,Q 2 )} , (16) 

again, P M refers to the nucleon momentum and Q 2 = —q 2 . The tensors multiplying the 
structure functions in eq. (0) should be identified with the Lorentz tensors e* in (f|). 
Contracting Wf$ with the longitudinal AjJ" and transverse Kjf projection operators [39 
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P • g 



9a 



r T t , 



P-g 



+ 2P- ga; 



^A + (g-5)g A P T e^ 



and choosing the pertinent polarization, yields the longitudinal component 



g L (x,Q 2 ) = gi(x,Q 2 



as well as the transverse combination 



g T (x, Q 2 ) = gi(x, Q 2 ) + g 2 {x, Q 2 ) . 



(17) 
(18) 

(19) 
(20) 



Also, b = —AM + 2P • qx — (q ■ S) 2 ^. In the Bjorken limit, which corresponds to 

the kinematical regime 



qo 



\q\ — M N x with \q\ — > oo , 



(21) 



the antisymmetric component of the hadronic tensor becomes |2S 

^ W k p sgn(fc ) 6 (k 2 ) dt 

x J d 3 xi J d 3 x 2 exp [—i(k + q) ■ (x\ — x 2 )\ 

x(N\ {ftxLtlQV^^.Ol + fKOlQYf*^,*)} \N) 



(22) 



where e Mpi/0 -7°"7 5 is the antisymmetric combination of 7 At 7 p 7,/. The matrix element between 
the nucleon states is to be taken in the space of the collective coordinates, A(t) (see eqs. fll2]) 
and (|15|)) as the object in curly brackets is an operator in this space. In deriving the expres- 
sion (|22"D the /ree correlation function for the intermediate quark fields has been assumed^ 
after applying Wick's theorem to the product of quark currents in eq. (||). [p9| . The use of 
the free correlation function is justified because in the Bjorken limit fl2~ID the intermediate 
quark fields carry very large momenta and are hence not sensitive to typical soliton mo- 
menta. This procedure reduces the commutator [J^xx, t), Jl(x 2 , 0)] of the quark currents 
in the definition (0) to objects which are merely bilinear in the quark fields. Consequently, 
in the Bjorken limit ( |2~T|) the momentum, k, of the intermediate quark state is highly off- 
shell and hence is not sensitive to momenta typical for the soliton configuration. Therefore, 



'Adopting a dressed correlation will cause corrections starting at order twist-4 in QCD |2 
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the use of the free correlation function is a good approximation in this kinematical regime. 



Accordingly, the intermediate quark states are taken to be massless, cf. eq. ([22|). 

Since the NJL model is originally defined in terms of quark degrees of freedom, quark 
bilinears as in eq. fl22|) can be computed from the functional 



(q(x)Q 2 q(y)) = J DqDq q(x)Q 2 q(y) exp J d 4 x' C 

j DqDq exp J d^x'd^y' 5 4 (x' — y')C 



i5a(x, y) 



+ a(x',y')q(x')Q 2 q(y')) (23) 

J / a(x,y)=0 



The introduction of the bilocal source a(x,y) facilitates the functional bosonization after 
which eq. (^3]) takes the form 

faA^T rA log - V)D + a(», y)&)) L Mrf • (24) 

The operator D is defined in eq. @. The correlation (q{x)Q 2 q(y)) depends on the angle 
between x and y. Since in general the functional (|23"D involves quark states of all angular 
momenta (I) a technical difficulty arises because this angular dependence has to be treated 
numerically. The major purpose of the present paper is to demonstrate that polarized 
structure functions can indeed be computed from a chiral soliton. With this in mind we will 
adopt the valence quark approximation where the quark configurations in (^) are restricted 
to the valence quark level. Accordingly the valence quark wave-function ([14]) is substituted 



into eq. ([22]). Then only quark orbital angular momenta up to I = 2 are relevant. From a 
physical point of view this approximation is justified for moderate constituent quark masses 
(m fa 400MeV) because in that parameter region the soliton properties are dominated 
by their valence quark contributions [2B|, [IB|. In particular this is the case for the axial 
properties of the nucleon. 

In the next step the polarized structure functions, gi(x,fi 2 ) and gT(x,fi 2 ), are extracted 
according to eqs. (|17]) and (18). In the remainder of this section we will omit explicit 



reference to the scale n 2 . We choose the frame such that the nucleon is polarized along the 
positive-z and positive-a; directions in the longitudinal and transverse cases, respectively. 
Note also that this implies the choice q = qz. When extracting the structure functions the 
integrals over the time coordinate in eq. ([22]) can readily be done yielding the conservation 
of energy for forward and backward moving intermediate quarks. Carrying out the integrals 
over k and k = \k\ gives for the structure functions 
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9l (x) = -N C ^(N, h\ J <m k k 2 \ $(p) (l-tx-k) 7 5 r^ v (p) 



k=qo+e v 



+#(-p) (l - a • fc) 7 5 r^ v (-p) t W iz) , (25) 

\ / k=q — e v 2 



9t(x) = g x (x) + g 2 (x) 

= -n c ^-(n, \x\ J dn k k 2 Ut( P ) (a • k) 7 5 r^(p; 



k=qo+£v 



+#(-p) (a ■ fe) 7 5 r> v (-p) . )\N, \x) , (26) 



where p = k + q and T = ^1 + ^D^Ti with .D^ = ~ tr ^TjA(t)rjA^j being the adjoint 
representation of the collective rotation cf. eq. (0). The second entry in the states labels 
the spin orientation. Nc appears as a multiplicative factor because the functional trace fl2"4l) 
includes the color trace as well. Furthermore the Fourier transform of the valence quark 
wave-function 

~ f d?x 

Vv(p) = J ipv{x) exp (ip ■ x) (27) 

has been introduced. Also, note that the wave-function ip v contains an implicit dependence 
on the collective coordinates through the angular velocity fl, cf. eq. (|14D. 

The dependence of the wave-function ip(±p) on the integration variable k is only implicit. 
In the Bjorken limit the integration variables may then be changed to El 



k 2 dVl k = pdp d& , p = \p\ , (28) 

where <3> denotes the azimuth-angle between q and p. The lower bound for the p-integral is 
adopted when k and q are anti-parallel; p± m = |Mjya;=Fe v | for k = — (q ± e v ), respectively. 
Since the wave-function ip(±p) acquires its dominant support for p < the integrand is 
different from zero only when q and k are anti-parallel. We may therefore take k = — z. This 
is nothing but the light-cone description for computing structure functions ||42| . Although 



expected, this result is non-trivial and will only come out in models which have a current 
operator which, as in QCD, is formally identical to the one of non-interacting quarks. The 
valence quark state possesses positive parity yielding ip(—p) = 7oV , (p)- With this we arrive 
at the expression for the isoscalar and isovector parts of the polarized structure function in 
the valence quark approximation, 

gl^ix) = -N c ^{N,h\ pdp d<$> 

187T 2 JM N \x T \ JO 

x^ v (p T ) (1 ± a 3 ) 7 5 Vv(p T )|iV, \z) (29) 
11 



Mat 1 f°° f 27r 

g{: ± \x) = -N C -^(N, -z\D 3i / pdp / d$ 

67T 2 JM w |x T | JO 



jv|a 

x^(p T ) Ti (i ± a3 ) T 5 vi v (p T )|iV,iz) , (30) 



^=0(x) = -N C ^-(N, -x\ / prfp / d$ 



x^(p>s7 5 Vv(p T )|iV,^> , (31) 



Mat 1 Z" 00 /' 27r 

gl=l( x ) = - N "( N , -a^ / Pdp / rf$ 

07T 2 JM N \x^\ JO 



M N \x T \ 

x^l(p T )T ia3l 5 Mp T )\N, -x) , (32) 

where x± = x ± e v /Mjv and cos(6p ) = M^x±/p. The complete structure functions are 
given by 



SiC*) = g{=°(x) + gi=\x) - (g{=°(x) - g[;_\x)) (33) 
9t(x) = g£*(x) + g^(x) - {g*?(x) - g^x)) . (34) 

Note also, that we have made explicit the isoscalar (J = 0) and isovector (1 = 1) parts. The 
wave-function implicitly depends on x because ip v (p±) = ip v (p,Q^,Q) where the polar- 
angle, 0^ , between p± and q is fixed for a given value of the Bjorken scaling variable x. 

Turning to the evaluation of the nucleon matrix elements defined above we first note that 
the Fourier transform of the wave-function is easily obtained because the angular parts are 
tensor spherical harmonics in both coordinate and momentum spaces. Hence, only the radial 
part requires numerical treatment. Performing straightforwardly the azimuthal integrations 
in eqs. (^) and ([3(]) reveals that the surviving isoscalar part of the longitudinal structure 
function, g[ =0 , is linear in the angular velocity, O. It is this part which is associated with 
the proton-spin puzzle. Using the standard quantization condition, Q = J / a 2 , where a 2 is 
the moment of inertia of the soliton and further noting that the z-direction is distinct, the 
required nucleon matrix elements are (N, \z\ J Z \N, \z) = ^. Thus, g[ =0 is identical for all 



nucleon states. Choosing a symmetric ordering , jl9| for the non-commuting operators, 
DiaJj \ {Di a , Jj} we find that the nucleon matrix elements associated with the cranking 
portion of the isovector piece, (N, ±|z| {D 3y , J x } \N, vanish. With this ordering we 

avoid the occurance of PCAC violating pieces in the axial current. The surviving terms stem 
solely from the classical part of the valence quark wave-function, \l/ v (x) in combination 
with the collective Wigner-D function, D 3z . Again singling out the ^-direction, the nucleon 



matrix elements become |T9] 
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(N,±z\D 3z \N,~z) = ~i 3 , (35) 

where i 3 = ±| is the nucleoli isospin. For the transverse structure function, the surviving 
piece of the isoscalar contribution is again linear in the angular velocities. The transversally 
polarized nucleon gives rise to the matrix elements, (N, ~x\J x \N, \x) = |. Again choosing 
symmetric ordering for terms arising from the cranking contribution, the nucleon matrix ele- 
ments (N, \x\ {D 3y , Jy] \N, \x) and (N, ~cc| {-D33, J y } \N, \x) vanish. As in the longitudinal 
case, there is a surviving isovector contribution stemming solely from the classical part of the 
valence quark wave-function, ^ v (x) in combination with the collective Wigner-D function, 
D 3x . Now singling out the ^-direction the relevant nucleon matrix elements become [B5|, 

(N 7 ~x\D 3x \N,~x} = ~i 3 . (36) 

Explicit expressions in terms of the valence quark wave functions (|33] and 0) for g["±(x), 
9i±{x)i 5 , 2.±°( a; ) an d Q 1 ^ 1 ^) are listed in the appendix A. 

Using the expressions given in the appendix A it is straightforward to verify the Bjorken 
sum rule |53j 

r? - r? = t dx {g{{x) - g n x {x)) = g A /Q , (37) 
Jo 

the Burkhardt-Cottingham sum rule |54| 

rf = f 1 dx g p 2 (x) = , (38) 

J 

as well as the axial singlet charge 

r; + r? = /' dx {&{ x ) + g?(x)) = 9 ° A , (39) 

Jo 

in this model calculation when the moment of inertia a 2 , as well as the axial charges g\ and 
gAi are confined to their dominating valence quark pieces. We have used 

AT r 

g A = ~ J rf 3 r^(r) 7 375r 3 ^v(r) (40) 
A/" r 

9a = -T d"ri,\{r) l3l ^{r) . (41) 

to verify the Bjorken Sum rule as well as the axial singlet charge. This serves as an analytic 
check on our treatment. Here a 2 , refers to the valence quark contribution to the moment 
of inertia, i.e. a 2 = (1/2) | (/^|t 3 |v) | 2 / (e^ — e v ). The restriction to the valence quark 
piece is required by consistency with the Adler sum rule in the calculation of the unpolarized 



structure functions in this approximation [47 
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V. NUMERICAL RESULTS 



In this section we display the results of the spin-polarized structure functions calculated 
from eqs. ( [A14jdAl~7|) for constituent quark masses of m — 400MeV and 450MeV. In 



addition to checking the above mentioned sum rules see eqs. (|37|) Cp9|) , we have numerically 
calculated the first moment of gf(x, /i 2 )| 

r? = f 1 dx g{{x) , (42) 
Jo 

and the Efremov-Leader-Teryaev (ELT) sum rule 

Tetl = f dx x {g{{x) + 2g»(x)) . (43) 

J 

We summarize the results for the sum rules in table 1. When comparing these results 
with the experimental data one observes two short-comings, which are already known from 
studies of the static properties in this model. First, the axial charge g& ~ 0.73 comes out 
too low as the experimental value is = 1.25. It has recently been speculated that a 



different ordering of the collective operators D ai Jj (cf. section 4) may fill the gap [55,46 
However, since such an ordering unfortunately gives rise to PCAC violating contributions 
to the axial current |48j and furthermore inconsistencies with G-parity may occur in the 



valence quark approximation |^9[ we will not pursue this issue any further at this time. 
Second, the predicted axial singlet charge g° A ~ 0.6 is approximately twice as large as 
the number extracted from experiment^ 0.27 ± 0.04 [17]]. This can be traced back to the 
valence quark approximation as there are direct and indirect contributions to g\ from both 
the polarized vacuum and the valence quark level. Before canonical quantization of the 
collective coordinates one finds a sum of valence and vacuum pieces 

9°A = 2(g° v + gL)^ = 4^¥- ( 44 ) 



3 Which in this case amounts to the Ellis-Jaffe sum rule [56] since we have omitted the strange 
degrees of freedom. A careful treatment of symmetry breaking effects indicates that the role of the 
strange quarks is less important than originally assumed [fL4,57|. 

4 Note that this analysis assumes SU(3) flavor symmetry, which, of course, is not manifest in our 
two flavor model. 
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Numerically the vacuum piece is negligible, i.e. g® ac / g® ~ 2%. Canonical quantization 
subsequently involves the moment of inertia a 2 = a 2 + a 2 ac , which also has valence and 
vacuum pieces. In this case, however, the vacuum part is not so small: a 2 ac /a 2 ~ 25%. 
Hence the full treatment of the polarized vacuum will drastically improve the agreement 
with the empirical value for g° A . On the other hand our model calculation nicely reproduces 
the Ellis-Jaffe sum rule since the empirical value is 0.136. Note that this comparison is 
legitimate since neither the derivation of this sum rule nor our model imply strange quarks. 
While the vanishing Burkhardt-Cottingham sum rule can be shown analytically in this 
model, the small value for the Efremov-Leader-Teryaev sum rule is a numerical prediction. 



Recently, it has been demonstrated [ 23 1 that that the ELT sum rule (f43|), which is derived 



within the parton model, neither vanishes in the Center of Mass bag model [23] nor is 
supported by the SLAC E143 data 0. This is also the case for our NJL-model calculation 
as can be seen from table I. 

In figure 1 we display the spin structure functions g{{x,^ 2 ) and gf^;^ 2 ) along with the 
twist -2 piece, g^ w ^ (x, /i 2 ) and twist-3 piece, ^(x^ 2 ). The actual value for /i 2 will be 
given in the proceeding section in the context of the evolution procedure. We observe that the 
structure functions g^x^jj 2 ) and g^ W<yP \x, /i 2 ) are well localized in the interval < x < 1, 
while for g\ about 0.3% of the first moment, Y\ = Jq dx g p (x,[i 2 ) comes from the region, 
x > 1. The polarized structure function g p (x, /x 2 ) exhibits a pronounced maximum at x ~ 0.3 
which is smeared out when the constituent quark mass increases. This can be understood 
as follows: In our chiral soliton model the constituent mass serves as a coupling constant 
of the quarks to the chiral field (see eqs. (|7|) and (|10"D). The valence quark becomes more 
strongly bound as the constituent quark mass increases. In this case the lower components 
of the valence quark wave-function increase and relativistic effects become more important 
resulting in a broadening of the maximum. With regard to the Burkhardt-Cottingham 
sum rule the polarized structure function g p (x,fi 2 ) possesses a node. Apparently this node 
appears at approximately the same value of the Bjorken variable x as the maximum of 
g p (x, /i 2 ). Note also that the distinct twist contributions to g p (%, yU 2 ) by construction diverge 
as In ( — > while their sum stays finite(see section 6 for details). 

As the results displayed in figure 1 are the central issue of our calculation it is of great 
interest to compare them with the available data. As for all effective low-energy models of 
the nucleon, the predicted results are at a lower scale Q 2 than the experimental data. In 
order to carry out a sensible comparison either the model results have to be evolved upward 
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or the QCD renormalization group equations have to be used to extract structure functions 
at a low-renormalization point. For the combination xgi(x) a parametrization of the em- 



pirical structure function is available at a low scale [0f]. In that study the experimental 
high Q 2 data are evolved to the low-renormalization point /i 2 , which is defined as the low- 
est Q 2 satisfying the positivity constraint between the polarized and unpolarized structure 
functions. In a next-to-leading order calculation those authors found /i 2 = 0.34GeV 2 |B2 



In figure 2 we compare our results for two different constituent quark masses with that 
parametrization. We observe that our predictions reproduce gross features like the position 
of the maximum. This agreement is the more pronounced the lower the constituent quark 
is, i.e. the agreement improves as the applicability of the valence quark approximation be- 
comes more justified. Unfortunately, such a parametrization is currently not available for the 
transverse structure function Qt{x) (or g 2 (x)). In order to nevertheless be able to compare 
our corresponding results with the (few) available data we will apply leading order evolution 
techniques to the structure functions calculated in the valence quark approximation to the 
NJL-soliton model. This will be subject of the following section. 



VI. PROJECTION AND EVOLUTION 



One notices that our baryon states are not momentum eigenstates causing the structure 
functions (see figures 1 and 2) not to vanish exactly for x > 1 although the contributions 
for x > 1 are very small. This short-coming is due to the localized field configuration and 
thus the nucleon not being a representation of the Poincare group which is common to the 
low-energy effective models. The most feasible procedure to cure this problem is to apply 



Jaffe's prescription | 33fl 



5 These authors also provide a low scale parametrization of quark distribution functions. However, 
these refer to the distributions of perturbatively interacting partons. Distributions for the NJL- 
model constituent quarks could in principle be extracted from eqs. (p?9|)-(32). It is important to 



stress that these distributions may not be compared to those of ref |32|] because the associated 
quarks fields are different in nature. 
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f(x) — f{x) = —f (-log(l - x)) (45) 

to project any structure function /(x) onto the interval [0,1]. In view of the kinematic 
regime of DIS this prescription, which was derived in a Lorentz invariant fashion within the 
1+1 dimensional bag model, is a reasonable approximation. It is important to note in the 
NJL model the unprojected nucleon wave-function (including the cranking piece^, see |14]) is 



anything but a product of Dirac-spinors. In this context, techniques such as Peierls-Yoccoz 
(which does not completely enforce proper support ||60|| , < x < 1 nor restore Lorentz 



invariance, see |6l[]) appear to be infeasible. Thus, given the manner in which the nucleon 
arises in chiral-soliton models Jaffe's projection technique is quite well suited. It is also 
important to note that, by construction, sum rules are not effected by this projection, i.e. 
/ °° dxf(x) = Jo dxf(x). Accordingly the sum-rules of the previous section remain intact. 

With regard to evolution of the spin-polarized structure functions applying the OPE 
analysis of Section 2, Jaffe and Ji brought to light that to leading order in 1/Q 2 , gi(x,Q 2 ) 
receives only a leading order twist-2 contribution, while g2(x,Q 2 ) possesses contributions 
from both twist-2 and twist-3 operators; the twist-3 portion coming from spin-dependent 
gluonic-quark correlations [4"0], [^T| (see also, |5D| and |)1J). In the impulse approximation 



40| , these leading contributions are given by 

lim [ l dxx n 9l (x,Q 2 ) = W °li > n = 0,2,4,..., (46) 

Q 2 -^oo JO 2 ^ 

g lto £ d x * g, (x , Q 2 ) = -fJL^ £ {0S> - Ot) , „ = 2, 4 (47) 

Note that there is no sum rule for the first moment, T 2 (Q 2 ) = Jq 1 dxg 2 (x, Q 2 ), J4(|. Sometime 
ago Wandzura and Wilczek [fj^] proposed that g 2 (x, Q 2 ) was given in terms of g±(x, Q 2 ), 



x,Q 2 ) = - 9l (x,Q 2 )+ f ^ 9l {y,Q 2 ) (48) 

Jx y 



which follows immediately from eqs. (f46j) and (^) by neglecting the twist-3 portion in the 



sum in (f47D. One may reformulate this argument to extract the twist-3 piece 

U*,Q 2 ) = 92(x,Q 2 ) - g^ w {x,Q 2 ) , (49) 



6 Which in fact yields the leading order to the Adler sum rule, F" p — F" p |47j] rather than being 
a correction. 
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since, 



j o dx x n g 2 {x,Q 2 ) = 2 ( w + i) ^°3,t , n = 2,4, .... (50) 

In the NJL model as in the bag-model there are no explicit gluon degrees of freedom, 
however, in both models twist-3 contributions to g 2 {x,fi 2 ) exist. In contrast to the bag 
model where the bag boundary simulates the quark-gluon and gluon-gluon correlations 
23| in the NJL model the gluon degrees of freedom, having been "integrated" out, leave 
correlations characterized by the four-point quark coupling Gnjl- This is the source of the 
twist-3 contribution to g 2 (x,fi 2 ), which is shown in figure 1. 

For g 1 (x, Q 2 ) and the twist-2 piece g 2 ww {x, Q 2 ) we apply the leading order (in olqcd(Q 2 )) 
Altarelli-Parisi equations |34| to evolve the structure functions from the model scale, fi 2 , to 
that of the experiment Q 2 , by iterating 

g(x, t + 6t) = g{x, t) + St^j^ , (51) 

where t = log (Q 2 / -^-qcd) ■ The explicit expression for the evolution differential equation is 
given by the convolution integral, 
dg(x,t) a(t) 



dt 2n 



■g(x,t) <g> P qq {x) 



"V Cr{f) f 1 «i Piqiy)g (* t ) (52) 



7T Jx y \y 

where the quantity P qq (z) = (-jzf?) represents the quark probability to emit a gluon such 

n' 2 — 1 

that the momentum of the quark is reduced by the fraction z. Cn(f) = ^- for n/-flavors, 

4-7T 

/31og(QVA 2 ) 



a QCD = flio g (Q2/A 2 ) and @ = ( n ~~ l n f)- Employing the "+" prescription @ yields 



dg(x,t) 2C R (f) [ ( x 2 \ 

" 1 x H h 2 log(l — x) ] g[x, t) 



dt 9 1 { \ 2 

,1 {l + y 2> 



+ L dy \- 



- 9 \ -,t ] -g{x,t) 
y \y 



y j yy \v j J J 

As discussed in section 2 the initial value for integrating the differential equation is given 
by the scale fi 2 at which the model is defined. It should be emphasized that this scale 
essentially is a new parameter of the model. For a given constituent quark mass we fit fi 2 
to maximize the agreement of the predictions with the experimental data on previously 
calculated unpolarized structure functions for (anti)neutrino-proton scattering: F% p — F 2 V - 
For the constituent quark mass m = 400MeV we have obtained [i 2 ~ 0.4GeV 2 . One certainly 



18 



wonders whether for such a low scale the restriction to first order in cxqcd is reliable. There 
are two answers. First, the studies in this section aim at showing that the required evolution 
indeed improves the agreement with the experimental data and, second, in the bag model 
it has recently been shown |62[ that a second order evolution just increases fi 2 without 



significantly changing the evolved data. In figure 3 we compare the unevolved, projected, 
structure function g\{x,jji 2 ) with the one evolved from /i 2 = 0.4GeV 2 to Q 2 = 3.0GeV 2 . 
Also the data from the E143-collaboration from SLAC [0 are given. Furthermore in figure 
3 we compare the projected, unevolved structure function ' (x,[i 2 ) as well as the one 

evolved to Q 2 = 5.0GeV 2 with the data from the recent E143-collaboration at SLAC 0. As 
expected we observe that the evolution pronounces the structure function at low x; thereby 
improving the agreement with the experimental data. This change towards small x is a 
general feature of the projection and evolution process and presumably not very sensitive 
to the prescription applied here. In particular, choosing an alternative projection technique 
may easily be compensated by an appropriate variation of the scale \j? . 

While the evolution of the structure function g 1 (x, Q 2 ) and the twist-2 piece g\ vw (x, Q 2 ) 
from /i 2 to Q 2 can be performed straightforwardly using the ordinary Altarelli-Parisi equa- 
tions this is not the case with the twist-3 piece g 2 (x, Q 2 ). As the twist -3 quark and quark- 
gluon operators mix the number of independent operators contributing to the twist-3 piece 
increases with n, where n refers to the n— moment foT| . We apply an approximation (see 
appendix B) suggested in [^] where it is demonstrated that in N c — > oo limit the quark 
operators of twist-3 decouple from the evolution equation for the quark-gluon operators of 
the same twist resulting in a unique evolution scheme. This scheme is particularly suited 
for the NJL-chiral soliton model, as the soliton picture for baryons is based on N c — > oo 
arguments^ 

In figure 4 we compare the projected unevolved structure function g 2 (x, jj?) evolved to 
Q 2 = 5.0GeV 2 using the scheme suggested in [j36|| . In addition we reconstruct g 2 (x, Q 2 ) at 
Q 2 = 3.0GeV 2 from g^ w ^ Q*) an d g 2 (x, Q 2 ) and compare it with the recent SLAC data 
for g 2 (x, Q 2 ). As is evident our model calculation of g 2 (x, Q 2 ), built up from its twist-2 
and twist-3 pieces, agrees reasonably well with the data although the experimental errors 
are quite large. 



7 This scheme has also employed by Song [23] in the Center of Mass bag model 



19 



VII. SUMMARY AND OUTLOOK 



In this paper we have presented the calculation of the polarized nucleon structure func- 
tions g\ (x, Q 2 ) and #2 {x, Q 2 ) within a model which is based on chiral symmetry and its 
spontaneous breaking. Specifically we have employed the NJL chiral soliton model which 
reasonably describes the static properties of the nucleon ||26|| , |pES| . In this model the cur- 
rent operator is formally identical to the one in an non-interacting relativistic quark model. 
While the quark fields become functionals of the chiral soliton upon bosonization, this fea- 
ture enables one calculate the hadronic tensor. From this hadronic tensor we have then 
extracted the polarized structure functions within the valence quark approximation. As the 
explicit occupation of the valence quark level yields the major contribution (about 90%) to 
the associated static quantities like the axial charge this presumably is a justified approxi- 
mation. When cranking corrections are included this share may be reduced depending on 
whether or not the full moment of inertia is substituted. 

It needs to be stressed that in contrast to e.g. bag models the nucleon wave-function 
arises as a collective excitation of a non-perturbative meson field configuration. In particular, 
the incorporation of chiral symmetry leads to the distinct feature that the pion field cannot 
be treated perturbatively. Because of the hedgehog structure of this field one starts with 
grand spin symmetric quark wave-functions rather than direct products of spatial- and 
isospinors as in the bag model. On top of these grand spin wave-functions one has to 
include cranking corrections to generate states with the correct nucleon quantum numbers. 
Not only are these corrections sizable but even more importantly one would not be able to 
make any prediction on the flavor singlet combination of the polarized structure functions 
without them. The structure functions obtained in this manner are, of course, characterized 
by the scale of the low-energy effective model. We have confirmed this issue by obtaining 
a reasonable agreement of the model predictions for the structure function g-y of the proton 



with the low-renormalization point parametrization of ref [32]. In general this scale of the 
effective model essentially represents an intrinsic parameter of a model. For the NJL-soliton 
model we have previously determined this parameter from the behavior of the unpolarized 
structure functions under the Altarelli-Parisi evolution [0. Applying the same procedure 
to the polarized structure functions calculated in the NJL model yields good agreement with 
the data extracted from experiment, although the error bars on g\ (x, Q 2 ) are still sizable. 
In particular, the good agreement at low x indicates that to some extend gluonic effects 
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are already incorporated in the model. This can be understood by noting that the quark 
fields, which enter our calculation, are constituent quarks. They differ from the current 
quarks by a mesonic cloud which contains gluonic components. Furthermore, the existence 
of gluonic effects in the model would not be astonishing because we had already observed 
from the non-vanishing twist-3 part of #2 {%, Q 2 ), which in the OPE is associated with the 
quark-gluon interaction, that the model contains the main features allocated to the gluons. 

There is a wide avenue for further studies in this model. Of course, one would like to 
incorporate the effects of the polarized vacuum, although one expects from the results on 
the static axial properties that their direct contributions are negligible. It may be more 
illuminating to include the strange quarks within the valence quark approximation. This 
extension of the model seems to be demanded by the analysis of the proton spin puzzle. 
Technically two changes will occur. First, the collective matrix elements will be more com- 
plicated than in eqs. (|3~5|) and (|36"D because the nucleon wave-functions will be distorted 



SU(3) .D-functions in the presence of flavor symmetry breaking [63|,16|. Furthermore the 
valence quark wave-function ( |T4] ) will contain an additional correction due to different non- 



strange and strange constituent quark masses When these corrections are included 

direct information will be obtained on the contributions of the strange quarks to polarized 
nucleon structure functions. In particular the previously developed generalization to three 
flavors pi allows one to consistently include the effects of flavor symmetry breaking. 
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APPENDIX A: SPIN POLARIZED STRUCTURE FUNCTIONS 



In this appendix we summarize the explicit expressions for the structure functions, eqs. 
([2D|-|3"2"D . The first step is to construct the eigenfunctions of the single particle Dirac Hamil- 
tonian in coordinate space. As the hedgehog ansatz (Q) connects coordinate space with 
isospace, these eigenfunctions are also eigenstates of the grand spin operator 
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G 



I H h - 

2 2 



(Al) 



which is the sum of the total spin J and the isospin r/2. The spin itself is decomposed 
into orbital angular momentum I and intrinsic spin er/2. Denoting by M the grand spin 
projection quantum number the tensor spherical harmonics which are associated with the 
grand spin may be written as yf- M (r). Note that these tensor spherical harmonics are 
two-component spinors in both spin and isospin spaces. Given a fixed profile function 
0(r) the numerical diagonalization of the Dirac Hamiltonian (|T0|) yields the radial functions 
#' + '%),/f etc. in (cf. ref |§) 



ip(o.+;i)(r)y 



G,M 
G,G+ 

G,M 
G+l,G-\ 



+ 



igf^ 2 \r)y t 



G+1,G 



f fl 




ig(°>-*)(r)y t 



J IX 




(A2) 
(A3) 
(A4) 



The second superscript (±) denotes the intrinsic parity, which also is a conserved quantum 
number|. Note that for the G = channel, which contains the classical contribution to the 
valence quark wave-function in eq. (TH 



Wv(r)y°i 



°.2 

0,0 , 



(A5) 



V Ur)y[ 

only the components with j = +1/2 are allowed. In addition to the classical piece ( |A5| ) 
the complete valence quark wave-function (|T¥) also contains the cranking correction, which 
dwells in the channel with G = 1 and negative intrinsic parity. 

The discretization (/^) is accomplished by choosing suitable boundary conditions at a ra- 
dial distance which is large compared to the soliton extension p5| , |64| . This calculation yields 



the energy eigenvalues e M , which enter the energy functional (11)- The soliton configuration 



is finally determined by self-consistently minimizing this energy functional. In ref. |66[ the 
numerical procedure is described in detail. 

We continue by making explicit the Fourier transform (|27|) of eq. (Tbl 



^ v (p) = ^ v (p) + Q M ^ M (p) . 



(A6) 



8 The total parity is given by the product of the intrinsic parity and (— ) G . 
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The leading order in N c valence quark contribution is just Fourier transform of 



/ ~a, (?) y°A (p) ^ 



^ 2 



(A7) 



/ 



and the cranking correction involves the Fourier transform of spinor with G = 1 and negative 
intrinsic parity 

\ 

(A8) 



( gf (p) (p) - ~gf (p) y'of (p) 



; 2 

(p) (p) - (p) (p 



2 

l.Af 



/ 



Here 3^ / (p) are the Fourier transforms of the tensor spherical harmonics associated with 



the grand spin operator (|AT|) . The Fourier transform for the radial functions is defined by 



dr r jiipr)^^) 



(A9) 



Here the index I of the Bessel function denotes the orbital angular momentum of the associ- 
ated tensor spherical harmonic. We have suppressed the grand spin index on the transforms 
of the radial wave functions for convenience. For purposes of notation we have introduced 



the quantity, which arises in analyzing the matrix elements (see eq. (H)) 



e v - e, 



V2 



V2 



(tlx - Q 2 ) - hi,o^o \ h 



G a ,l 



(A10) 



where 



e v - 



dr r 2 {g v (r) gf (r) + / v (r) /f (r) } 



(All) 



Defining the following combinations, 



/ (t) (p) = (p) 



9 { Hp) =Q»~gf(p) , « = i,2 

the isoscalar (vector) contributions to the spin polarized structure functions, eqs. 
read 



(A12) 
(A13) 
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POO 



5M N r 

' Nc ^~ / P d P 

Son Jm n \x t \ 

< {g v (p) (p) 1 " 3c ^ 2(9p±) - g v (p) g (2) (p) \ 

^ T ~9AP) / (2) (P) 



=f/1 



cos(e^ 



^P)/ (1) (p)^f 

(P)^ (1) (P)^®T/ V (P)^ 2) (P) 

-a (p) (p) 1 + C ^ (9p±) + / v (p) ; (2) (p) ■ 



• 2cos 2 (9 J : 



(A14) 



367T JAfvlajil 



367T Afjv|a; T | 

(p) 2 ± 2~g v (p) / v (p) cos(9j ) - / v (p) 2 (l - 2cos 2 (6j)) } , (A15) 



x 



#£ ± °(x,/z 2 ) 



-N, 



5M, 



A' 



367T JM N \x 



pdp 



JVF=F 



r (1) 3cos 2 (e±)-i _ (2) i 

x j<7v G»)# (1) (p) — ^ # v G»)# (2) (p)g 



T2v (P) (?) 



cos2(0?) - 3 + / v (p)^(p) cos2(e - 



4^ 



(A16) 



-AT, 



Mat r°° 



c: 



pdp [g v (p) 2 - f v (p) 2 cos 2 (9j)} . 



, r- • ..•/> J • .. > •/ " ' I ' (A 17) 

367T iMjv|x T | 1 

which we evaluate numerically. Note that in case of the neutron the signs of the isovector 
pieces have to be reversed. Note that the angle 0^ is related to the integration variable p 



via 



cos6^ = - \M N x =F e v |. 



(A18) 



In ref |T7| structure function g\ was calculated omitting the cranking corrections. In the 
special case of the isovector component these corrections drop out and we formally confirm 
the result displayed in eq (B.6) of ref [|67 . 



APPENDIX B: EVOLUTION OF g 2 (x, \i 2 



In the N c —>■ oo limit it has been shown that one can evolve the moments of g 2 {x, /J, 2 ), 
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b 



MAQ 2 )=\^^r\ MA i?) (Bl) 



from the scale, fi 2 to Q 2 , where the anomalous dimensions are 

7,-i = 2iV c + ^ + 7e - ^) , (B2) 

with ip(x) = (d/dx) log T(x) and b = (11N C — Inf) /3. N c and n/ are the number of colors 
and flavors respectively. Given the moments of g 2 (x, fi 2 ) 

Mj (/i 2 ) = jT 1 dx x j -% (x,fj, 2 ) , (B3) 

and expressing ~g 2 (x, /i 2 ) in terms of the log (a;) and a power series in x, 

oo 

g 2 (x, fi 2 ) = ai (/i 2 ) log(x) + J2 a n (V ) x n (B4) 

n=0 

one can alternatively express the moments Mj in terms of the coefficients a n 

Mj (/i 2 ) = A jn a n (/i 2 ) . (B5) 

We calculate the moments, Mj (/j?) from ( P33|) and evolve them according to (|B 1|) . Finally, 
inverting the matrix Aj n we obtain the evolved coefficients, a n (Q 2 ) which in turn yields 
~g 2 (x, Q 2 ) (see figure 4). 
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TABLES 

TABLE I. Sum rules calculated from eqs. (38-0S|) as functions of the constituent quark mass 
m in the NJL chiraEsoliton model. 



m (MeV) 400 450 

Burkardt-Cottingham: 

Bjorken: T{ - BJ = g A /6 0.121 0.118 

Ellis-Jaffe: T{ 0.149 0.139 

ELT: Belt 1.38x10~ 2 7.65 xl0~ 3 

Axial Singlet Charge: B^ + B™ = g° A 0.638 0.579 
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Figure Captions 



Figure 1: The valence quark approximation of the polarized proton structure functions 
as a function of Bjorken-x. Left panel: g\ (x, /i 2 ) for two constituent quark masses m. 
Right panel: g\ (x, fi 2 ) (solid line), g^ w ^ (x, (long-dashed line) and twist three portion, 
Tf 2 {x,fi 2 ) (dashed line). In this case we have used m = 400MeV. 

Figure 2: The valence quark approximation to the nucleon structure function xg\(x) in the 
NJL-soliton model compared to the low-renormalization point result of ref [3^ . 



Figure 3: The projection and evolution of the spin-polarized structure functions as a 
function of Bjorken-x. For the constituent quark mass we choose m = 400MeV. Left 
panel: gf(x,Q 2 ), unprojected (long-dashed line), projected (dashed line) and evolved from 
\i = 0.4GeV 2 to Q 2 = 3.0GeV 2 (solid line). Data are from |7j. Right panel: g? W{p) (x,Q 2 ), 
unprojected (long-dashed line), projected (dashed line) and evolved from fi = 0.4GeV 2 to 
Q 2 = 5.0GeV 2 (solid line). Data are from and || , where the diamonds, circles and trian- 
gles correspond to the 4.5° E143, 7.0° E143 and SMC kinematics respectively. Overlapping 
data have been shifted slightly in x. The statistical error are displayed. 

Figure 4: The evolution of ~g\ (x, Q 2 ) (projected) from jj, = 0.4GeV 2 (long-dashed line) 
to Q 2 = 5.0GeV 2 (solid line). In addition we display the corresponding evolution for 
g Y w(p) (x, Q 2 ) (projected). Right panel, g2{x,Q 2 ) = g^ W(j ^ (x, Q 2 ) +g2{x,Q 2 ) evolved 
from /i 2 = 0.4GeV 2 to Q 2 = 5.0GeV 2 . Data and statistical errors for g\ (x, Q 2 ) are displayed 
from [Q], where the diamonds and circles correspond to the 4.5° E143, 7.0° E143 kinematics 
respectively. Overlapping data have been slightly shifted in x. 
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